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Abstract 

Rough surfaces submerged in a liquid can remain almost dry if the liquid does not fully wet the roughness 
and gases are sustained in roughness grooves. Such partially dry surfaces can help reduce drag or enhance 
boiling. Gases sustained in roughness grooves would be composed of air and the vapor phase of the liquid 
itself. The thermodynamics of sustaining vapor was considered in a prior work [Patankar, Soft Matter^ 
2010, 6, 1613]. Here, the thermodynamics of sustaining gases (e.g. air) is considered. Governing equations 
are presented along with a solution methodology to determine a critical condition to sustain gases. The 
critical roughness scale to sustain gases is estimated for different degrees of saturation of gases dissolved 
in the liquid. It is shown that roughness spacings of less than a micron are essential to sustain gases on 
surfaces submerged in water at atmospheric pressure. This is consistent with prior empirical data. 


1 Introduction 


Rough surfaces that exhibit non-wetting properties under submerged conditions are desirable for many 
applications like drag reduction, 0i boiling, 1^ among others. Here, non-wetting behavior is dehned as 
the one where the liquid (e.g. water), into which the rough surface is submerged, does not fully wet the 
surface. For this to happen, it is essential to sustain gases in roughness grooves of the surface. However, 
sustaining gases in roughness grooves over long time periods has been challenging. Typically, gases are 
found to deplete after 2-3 days from rough surfaces with tens of micron scale features, whereas 

submicron scale roughness is found to sustain gases for more than 120 days. Increased pressure is 
found to adversely affect the ability to sustain gas. An understanding of the governing equations is 
essential to design surfaces that can sustain gases over an extended period. [8-10 
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In general, gases sustained in roughness grooves would be a mixture of air and the vapor phase of the 


submerging liquid. Thermodynamics for sustaining vapor has been studied before. 11-14 In this work, 
thermodynamics of sustaining air is considered. Specifically, the goal is to present an analysis of the 
equilibrium state where a gas like air prefers to remain in roughness grooves under submerged conditions. 
In this analysis it is assumed that vapor produced by the submerging liquid is not present in roughness 


grooves. However, generalization is possible by using the concept of partial pressure. 15 The presence 
of the vapor would make the surface more non-wetting. 

In the next section, background on the thermodynamics of gas dissolution is presented. Following 
that, the theory is applied to an example problem of sustaining air in cylindrical pores. 


2 Thermodynamics of gas dissolution 


2.1 Henry’s law 


Consider a gas in chemical equilibrium with a liquid bath into which this same gas is dissolved (Fig. [^. 
Consider the liquid solvent to be non-volatile, i.e., it does not evaporate. This implies that there will be 
only one type of gas present above the liquid bath. It is generally observed that the greater the pressure of 
the gas, the greater will be the dissolved mole fraction of the gas in the liquid. The equilibrium pressure 
(of the gas outside the liquid) versus the mole fraction (of the gas dissolved in the liquid) relationship is 
given by Henry’s law: 



( 1 ) 


where Pg is the pressure of the gas outside the liquid, Xg^\ is the mole fraction of the gas dissolved in the 
liquid, and is Henry’s constant. Eqn. j^is valid for ideal solutions. Real solutions have non-idealities 
due to which the pressure to mole fraction ratio is not constant in general. 


2.2 Chemical equilibrium 

An expression for the pressure-mole fraction relationship at equilibrium is obtained by equating the 
chemical potential of the gas outside the liquid to the chemical potential of the gas dissolved in the 
liquid. 
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Figure 1. Chemical equilibrium between a gas (denoted by circles) dissolved in a liquid solvent and 
the same gas outside the solvent. The liquid solvent is assumed to be non-volatile. 


The chemical potential /ig of the gas outside the liquid is given by 


15 


16 


Mg = Mg[7"] + RT ln[V’gPg], 


( 2 ) 


where [•] implies “fundion-of,” is the standard potential of the gas which depends only on temperature 
T, R is the universal gas constant, ipg is the fugacity coefficient of the gas that accounts for non-idealities, 
and Pg is the pressure of the gas. In general, ijjg is not constant. 

The chemical potential fig^i of the dissolved gas is given by 


15 


16 


Mg.i = Mg.i[Pi,r] -F RT ln[7g,ia;g,i]. 


( 3 ) 


where 7 gj is the activity coefficient, Xg^\ is the mole fraction, and is the standard potential of 

the dissolved gas in the liquid that has pressure pi and temperature T. In general, 7 g_i is not constant. 

At equilibrium the chemical potentials of the gas phase and the dissolved gas must be equal. It follows 
from Eqns. and that 


a^g.i i’s 


[T] 


( 4 ) 


Eqn. I^is a general expression where the pressure Pg of the gas phase and the pressure p\ of the liquid need 
not be equal. For the case depicted in Fig. pg and p\ are indeed equal since liquid solvent is assumed 
to be non-volatile. However, consider another situation where there are multiple gas species above the 
liquid in Fig. [l] In this case, Pg would be the partial pressure of the gas of interest and pi would be equal 
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to the total pressure exerted by all gases on the liquid. 

If the gas is assumed to be ideal, then ipg = 1. Additionally, if the solution is dilute, which is often 
the case with gases dissolved in liquids, then 7 g_i —1. Using this in Eqn.|^ 


^ = H[pI^ 


,T]e 


Mg j[pf^ 


,T] 


( 5 ) 


where implies that pi = 1 atmosphere (atm.), which is chosen to be the reference pressure. Addi¬ 

tionally, 




( 6 ) 


It follows from Eqn. that Pg/xgj = when pi = 1 atm. Hence, is Henry’s 

constant (see Eqn. at 1 atm. pressure of the liquid. It is seen from Eqn. that the pressure to mole 
fraction ratio Pg/xg^i is dependent on the temperature, as well as the pressure of the liquid. 16 Eqn. 


can be simplified further to obtain 15 


^ = H[Pi 


latm^ ^ 


( 7 ) 


where Ugj is the partial molar volume of the gas dissolved in the liquid. In general, is dependent 
on a;g 1 , T, and pi. However, the variation of Ugj is not significant even upto liquid pressures as high as 


1000 atm. 16 


2.3 Isobaric equilibrium 

For the case depicted in Fig. ^ Pg = Pi = Pe, where Pe is the equilibrium pressure. Since the gas and 
liquid pressures are equal, this will be regarded as an isobaric equilibrium. For isobaric equilibrium, 
Eqn. [^becomes 


2^g.i = 


-u, 

Pee RT 






( 8 ) 


Eqn.j^gives the pressure-mole fraction curve for isobaric equilibrium. This curve, plotted in Fig.[^ will be 
called the isobaric equilibrium curve. Fig. [^specifically plots the case of oxygen (gas) dissolved in water 
(liquid), where = 25°C] = 42,590 atm and V^.i = 32 ml/mol. 16 Eqn. [^ whereassumptions 

for an ideal gas and dilute solution were made, will be less accurate at high mole fractions and high 
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pressures. Henry’s law is also plotted in Fig. with = 25°C] = 42,590 atm. It is seen 

that Henry’s law is a reasonable approximation to the isobaric equilibrium curve at low mole fractions and 
pressures. Eqn. [^implies that the isobaric equilibrium pressure Pe depends on T and Xgj, i.e., Pe[xg^\,T]. 



Figure 2. The isobaric equilibrium curve (Eqn. ^ compared to the equilibrium relation based on 


Henry’s law (Eqn. 1 
= 42, 590 atm, ancP 


, The gas is oxygen and the solvent is water, where = H[p\°‘^'^,T = 25°C] 
gj = 32 ml/mol. 


16 


2.4 Non-isobaric equilibrium 


Equilibrium in a general case where Pg ^ pi can be written in terms of the isobaric equilibrium pressure 
Pe[xg^\,T] as follows 


Ps _ g^Cpi-Pobg.i.T])^ 

Pe[Xg,hT] 


(9) 


Eqn. 1^ is obtained by using Eqns. [^and To understand Eqn. [^consider the following example. 

Imagine the same configuration as in Fig. but with a solution pressurized by a mixture of gases 
instead of only one gas. This configuration is one example of a non-isobaric equilibrium (a second example 
will be seen in Section 3, where submerged cylindrical pores are considered). The pressure, pi, of the 
liquid solvent will be the same as the total pressure exerted by all gases on the liquid. Let Xgj = xa 
be the mole fraction of one particular gas dissolved in the liquid at pressure pi — denoted by point A in 
Fig. § At equilibrium, what would be the partial pressure Pg of this particular gas in the gaseous mixture 
above the solution? The answer to this question is provided by Eqn. and denoted by point B in Eig. 
Since pi ^ Pg in this case, the chemical equilibrium between the dissolved gas (point A in Eig. and the 
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gas phase (point B in Fig. [^, will be regarded as a non-isobaric equilibrium. 

Given p\ and xx (point A in Fig. [^, the gas partial pressure Pg (point B in Fig.[^ can alternately be 
obtained by using Eqn. Given xa (and T ), the equilibrium pressure Pe [xa] for an isobaric configuration 
(Fig.[T]) is given by Eqn.[^ This is denoted by point G in Fig.[^ Using known values of pi (point A) and 
Pe[a:^A] (point G), the partial pressure of the gas Pg (point B) can be calculated according to Eqn. 

It is noted that Pg yf Pe{xp^. If p\ — Pe is not large, then pg (point B) may be approximated by Pe[a:A] 
(point G), which in turn may be approximated by pn at point D by using Henry’s law (pn = H^^xa] 
Fig.§. 



Figure 3. A schematic diagram to understand non-isobaric equilibrium between a gas dissolved in a 
liquid and the same gas adjacent to that liquid. 


2.5 Analogy between solubility and phase equilibrium calculations 

Gonsider Fig. to understand the analogy between solubility and phase equilibrium calculations. The 
pressure-mole fraction curve for isobaric equilibrium is analogous to the pressure-temperature co-existence 
(binodal) curve in liquid-vapor phase equilibrium. Eqn. is analogous to the Clasius-Glapeyron equation 


(see Eqns. 1-3 in Patankar, 2010 11 for example). Thus, the isobaric equilibrium pressure Pe in solubility 
calculations is analogous to the saturation pressure Psat in phase equilibrium calculations. 

A gas dissolving into a liquid in case of solubility is analogous to the condensation of a vapor to 


its liquid in case of phase change. Similarly, the release of a dissolved gas from a liquid (solubility) is 
analogous to the formation of vapor during boiling (phase change). 
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In case of solubility, dissolution of gas is favored above the isobaric equilibrium curve where the liquid 
is undersaturated, while release of gas from liquid is favored below the isobaric equilibrium curve where 
the liquid is supersaturated. Analogously, in case of phase change, condensation is favored above the 
co-existence curve where the liquid is stable, while boiling is favored below the co-existence curve where 
the liquid is metastable. 

Equating chemical potentials of the gas phase and the dissolved gas (Eqn. in solubility calculations 
is analogous to equating chemical potentials of the liquid and vapor phases (see, for example, Eqn. 4 in 


Patankar, 2010 11 ) in phase equilibrium calculations. 

Fig.j^is a visual aid to organize the calculation process. Specifically, point A with coordinates (pi, xa) 
represents the dissolved gas at mole fraction xa in a liquid at pressure pi. Point B with coordinates (pg, xa) 
represents the gas phase at pressure Pg that is in chemical equilibrium with the same gas dissolved at 
mole fraction xa in the liquid. As long as these interpretations are recognized. Fig. and the analogy 
with phase equilibrium, provides a useful framework to conceptualize the calculations. This will be made 
evident in subsequent sections. 


3 Sustaining gases on rough surfaces 

3.1 Calculation methodology 

Consider a rough surface with cylindrical pores that is completely submerged in water. Is it possible 
to sustain a gas in the pores? This question will be analyzed in this section. In this analysis it will be 
assumed that the length scales are below the capillary length scale of the liquid 2.72 mm for water) 
so that the effect of gravity is not dominant. 

Consider the cross-section of a typical cylindrical pore on a rough surface as shown in Fig. Let 
there be a liquid with a dissolved gas above the pore and the same gas inside the pore. For the gas 
to be sustained in the pore, there should be thermodynamic equilibrium which constitutes chemical, 
mechanical, and thermal equilibria. The equilibrium conditions can be derived from energy minimization 
similar to prior work, [11| which will not be presented here. The temperature will be assumed uniform 
everywhere so that thermal equilibrium is ensured. Chemical equilibrium implies that the gas in the pore 
should have the same chemical potential as the gas dissolved in the liquid above the pore. For simplicity, 
the liquid will be assumed to be non-volatile. Otherwise chemical equilibrium between liquid and vapor 
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Figure 4. Cross-section of a cylindrical pore on a surface submerged in water. Circles denote gas 
molecules either dissolved in the liquid or inside the pores. 


phases of the solvent should be considered, fl Mechanical equilibrium implies that the pressure difference 
between the liquid and the gas in the pore is balanced by the surface tension of the curved interface and 


that the curved interface remains pinned at the top of the pore (Fig. ffl. 11 


Consider Case 1, where the liquid pressure above the pore is p\. Let there be isobaric equilibrium 
between the gas dissolved in the liquid and the gas inside the pore. For isobaric equilibrium, Eqn. 
gives the mole fraction of the dissolved gas in the liquid (xg^i = Xsatipi]) and the gas pressure inside the 


pore would be equal to the liquid pressure, i.e., Pg = pi (Section 2.3). This equilibrium is denoted by 
point E in Fig. In this case, the liquid is considered to be saturated with the gas. The dissolved gas 
and the gas inside the pore would be in chemical equilibrium. Since the liquid and gas pressures are the 
same, mechanical equilibrium at the interface requires the liquid-gas interface to be flat. If the material 
of the surface is hydrophobic (i.e. contact angle > 90°) then this flat interface will remain pinned at the 


corner at the top of the pore 17 irrespective of the pore diameter as long as gravitational effect is not 
dominant. In this case, the gas will be sustained in the pore. 

Now consider Case 2, where the liquid pressure is the same as in Case 1, i.e. equal to pi, but the 
mole fraction of the dissolved gas is Xgj = xa < a;sat[pi]- In this case the liquid is undersaturated with 
dissolved gas (point A in Fig. [^. Let the degree of saturation </) of the dissolved gas denoted by point A 
in Fig.[^be quantified hy (j) = Xgj/xsat = a;A/a;sat[pi]- By imposing chemical equilibrium, the pressure Pg 
of the gas inside the pore can be calculated according to Eqn. [^or Eqn. This is denoted by point B in 
Fig. i The dissolved gas (point A) will be in chemical equilibrium with the gas in the pore (point B). 
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This is a case of non-isobaric equilibrium. Mechanical equilibrium requires that the liquid-gas interface 
will be curved and the radius of curvature i?co (Fig. of this interface is given by the Young-Laplace 
equation: 


Rco = 


2a\g 


( 10 ) 


(Pi - Ps) ’ 

where dig is the liquid-gas interfacial tension. Mechanical equilibrium also requires the liquid-gas interface 
to remain pinned at the top of the pore in order to sustain gas inside the pore. For this to happen the 


radius R of the cylindrical pore must satisfy the following pinning condition 11 


R < i?cr, where 

Rcr — .^coCOS^e- (^^) 


9e is the equilibrium material contact angle of the surface. Thus, the pore radius should be smaller than 
the critical radius i?cr- For Case 1, where the liquid is saturated with the gas, the critical pore radius is 
infinitely large (under the assumption that gravity is not important). For Case 2, there is a finite value of 
the critical pore radius. The critical radius will depend on the type of gas dissolved and other conditions 
such as liquid pressure and the degree of saturation of the dissolved gas, among others. It is, however, 
possible to simplify the governing equations and obtain estimates for critical pore sizes. This will be 
discussed in the next section. 

It can be verified from the governing equations that as the degree of saturation decreases, i.e., as 
(f> decreases, the pressure difference, p\ — Pg, between the liquid and the gas in the pore will increase 
for chemical equilibrium to be maintained. This implies smaller radius of curvature Rco of the liquid- 
gas interface. Consequently, the critical radius of the cylindrical pore would be smaller for increasingly 
undersaturated liquid. 

3.2 Approximate calculation of the critical pore size 

In this section, approximate calculations for Case 2 of the previous section are presented. The configu¬ 
ration is as in Fig. and calculations will be done with reference to Fig|^ 

Fig. 1^ shows the isobaric equilibrium curve and its comparison with Henry’s law for oxygen. Fig. 
shows similar curves for nitrogen and carbon dioxide - two other gases that are also present in air. Carbon 
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dioxide dissolves in much greater quantity in water compared to oxygen and nitrogen. It is seen from 
Figs. and that for all these gases Henry’s law is close to the isobaric curve for pressures as high 
as 100 atm. Hence, in the following calculations Henry’s law will be used to approximate the isobaric 
equilibrium curve for all gases - oxygen, nitrogen, and carbon dioxide. In the interest of simplicity, in 
the following calculations, each of these gases will be considered independently and not as a mixture. 




Figure 5. Isobaric equilibrium curves (Eqn. compared to the equilibrium relation based on Henry’s 
law (Eqn. [T|) for nitrogen (top) and carbon dioxide (bottom). = 90,770 atm and Egj = 33.3 ml/mol 


for nitrogen, 16 and 77° [ = 1,630 atm and Vg^i = 37.6 ml/mol for carbon dioxide. 
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As noted at the end of Section [2^ the gas pressure in the pore, which is denoted by point B in Fig.[^ 
will be approximated by point D in Fig. The error is typically not more than 15%. For example, in 
the case of oxygen, it can be verified that the error in calculating oxygen pressure in the pore (difference 
between pg at point B and pn at point D in Fig. is around 14% when it is in equilibrium with oxygen 
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that is dissolved in water at 100 atm. (point A in Fig.[^. While, in the current application the error may 


be disregarded, this difference cannot be ignored for deep sea organisms and deep diving mammals. 16 


Using approximations above, the following equations are to be noted. First, the degree of saturation 
(j) of the gas dissolved in the liquid (point A in Fig. above the pore (Fig. can be approximated by 


0 


XA 

xh’ 


( 12 ) 


where a;H is the saturation mole fraction based on Henry’s law (point F in Fig. [^. a;H is used to 
approximate the saturation mole fraction a;sat[pi] at point E in Fig. Note that the error between a;H 
and Xsatbi] is exaggerated in the schematic in Fig.[^ 

Second, xh is determined using Henry’s law (point F in Fig. as follows 


XH = 


P\ 


H: 


(13) 


g.i 


Third, the pressure of the gas in the pore can be approximated by (refer to Fig. 


Pg « Pe[xA] « PH = XA- 
Eqns. [T2p4| lead to the following relation 


(14) 


Xa 

xh 


Ps 
P\' 


(15) 


Use of Eqn. ]_5 in Eqns. W and El leads to the following equation for the critical pore radius 


Rcr «-^^^^^COS^e- (16) 

Pl(l - 0) 

The cylindrical pore radius should be less than the critical value obtained from Eqn. to sustain gas 
within a cylindrical pore. Note that, with the approximations used, the critical radius does not depend 
on the type of gas being considered once the degree of saturation (j) is known. 

Eig.§ shows the critical pore radius for different degrees of saturation of gases dissolved in water at 
Pi = 1 atm. The surface tension of the water-gas interface is a\g = 72 mN/m at room temperature and 
0 = 110° is considered (typical for hydrophobic chemical coatings on rough surfaces). Fig. shows that 
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Figure 6. Critical pore radius to sustain gas in cylindrical pores on a surface that is submerged in 
water. Pore radii should be below the critical radius to sustain gas. (p is the degree of saturation of the 
gas dissolved in liquid water. Liquid water is assumed to be at 1 atm., cjig = 72 mN/m, and 0 = 110°. 


the smallest critical pore radius at 1 atm. water pressure is 480 nm corresponding to ^ = 0, i.e. no 
dissolved gases in water. In other words pore diameters of less than a micron are necessary to prevent the 
liquid from wetting the pores irrespective of the degree of saturation of dissolved gases. Similar analysis 


can be extended to pillar or other geometries without fundamental difficulty. 12 Calculations predict 
that pillar spacing of less than a micron would be required to prevent wetting of rough surfaces submerged 
in water. This is consistent with experimental results where it is observed that surfaces remain dry under 
water when roughness length scales are below micron scale, whereas surfaces with tens of micron scale 
roughness typically get wet in 3-4 days for liquid pressures around 1 atm. 

It also follows from Eqn. that the critical radius decreases with increasing liquid pressure for the 
same degree of saturation. For example, the critical pore radius would reduce from 480 nm to 48 nm if 
the water pressure were increased from 1 atm. to 10 atm. 


4 Conclusions 

Thermodynamics of sustaining gases (e.g. air) on rough surfaces submerged in water is considered. 
Governing equations and a solution methodology to determine critical condition to sustain gases is pre¬ 
sented. It is shown that roughness spacings of less than micron are essential to sustain gases on surfaces 
submerged in water at atmospheric pressure. This is consistent with prior empirical data. 
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In the analysis presented here, the liquid, into which the rough surfaces would be submerged, was 
considered to be either saturated or undersaturated with dissolved gases. Extension to supersaturated 
cases can be done without difficulty. In that case the liquid pressure would be lower than the gas 
pressure at equilibrium. Consequently, sustaining gases would be possible even with hydrophilic surfaces. 
However, if the liquid is undersaturated with gas then hydrophilic surfaces submerged in it would not be 


able to sustain gases inside cylindrical pores - reentrant roughness geometry would be necessary. 19, 20| 
In general superhydrophilic surfaces would promote wetting of submerged surfaces. 
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